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ABSTRACT 

We examine the origin of clustercentric gradients in the star formation rates and colors of rich cluster 
galaxies within the context of a simple model where clusters are built through the ongoing accretion of 
field galaxies. The model assumes that after galaxies enter the cluster their star formation rates decline 
on a timescale of a few Gyrs, the typical gas consumption timescale of disk galaxies in the field. Such 
behaviour might be expected if tides and ram pressure strip ofl[ the gaseous envelopes that normally fuel 
star formation in spirals over a Hubble time. Combining these timescales with mass accretion histories 
derived from N-body simulations of cluster formation in a ACDM universe, we reproduce the systematic 
differences observed in the color distribution of cluster and field galaxies, as well as the strong suppression 
of star formation in cluster galaxies and its dependence on clustercentric radius. The simulations also 
indicate that a significant fraction of galaxies beyond the virial radius of the cluster may have been 
within the main body of the cluster in the past, a result that explains naturally why star formation 
in the outskirts of clusters (and as far out as two virial radii) is systematically suppressed relative to 
the field. The agreement with the data beyond the cluster virial radius is also improved if we assume 
that stripping happens within lower mass systems, before the galaxy is accreted into the main body of 
the cluster. We conclude that the star formation rates of cluster galaxies depend primarily on the time 
elapsed since their accretion onto massive virialized systems, and that the cessation of star formation 
may have taken place gradually over a few Gyrs. 



1. INTRODUCTION 



1999). Tides, either by the main cluster potential (e.g. 



Extensive observational work has established that the 
star formation properties of galaxies in rich cl usters differ 
significantly from those o f field galaxies (e.g 



Ostcrbrock 



1960; Dressier ct al. 1985; Balogh et al. 1997 



Koopman: 



i 



and Kenney 1998 ). In the field, galaxies form stars at rates 
several times higher than systems of similar luminosity in 
the cores of clusters. This is partly a result of the well- 
known morphology-density relation, since ellipticals and 
SO galaxies are more ab undant in clusters (Dressier 1980 



|Byrd and Vahonen 1990| ; [Fujita 19981) or by fiy-by en 
counters with othe r cluster galaxies ( "galaxy harassment" , 
Moore et al. 1996), have also been proposed to explain the 
observations. 

Although the above processes are all plausible, the ac- 
tual changes in star formation rate induced by them re- 
main a matter of debate. For example, some authors 
have suggested that the cluster environment triggers in- 
tense bursts of star formation t hat rapidly consume the 



Whitmorc ct al. 1993), but there is evidence that even 



later type galaxies in clusters form stars at lo wer r ates 
then in the field. For example, Balogh et al ( 199£ ) re- 
port that field galaxies of given bulge-to-disk ratio and 
luminosity have, on average, much larger [OH] equivalent 
widths than their counterparts in rich clusters, suggest- 
ing that the cluster environment somehow curbs the star 
formation rates of all galaxies, regardless of morphology. 

Various physical mechanisms have been proposed to ex- 
plain this and other systematic differences between the 
field and cluster galaxy populations. Ram pressure strip- 
ping by the intracluster medium, for example, has been 
suggested as a means of removing the gaseous component 
of disk galaxies and of dramatically altering their mor 



gas of an infalling cluster galaxy (Dressier and Gunn 1983 
Bargcr et al. 199^ ; Poggianti et al. 1999 ), while others 



have favored a scenario where star formation is truncated 
in a galaxy almost immediately after it is accreted into 
the cluster (Abraham ct al. 1996; Newberry et al. 1990 



Morris et al. 1998|; [Jones et al. 2000|y 



The virtues of these scenarios can in principle be tested 
observationally, since they are expected to produce a popu- 
lation of galaxies with unusually strong Balmer absorption 
lines in which star formation has been recently terminate d 
(e.g.. Dressier and Gunn 1983: Couch and Sharpies 1987). 



Recently, Balogh et al. (1999) have applied this idea to 
galaxies in the C NOCl survey of X-ray luminous clusters 
( Yec et al. 1996 ). From the paucity of galaxies with strong 



phology and subsequent star forming history (e.g.. Got! 



and Gunn 1972; Fujita and Nagashima 1999; Abadi et al 



H(5 absorption in their spectra they conclude that the de- 
cline of star formation in these clusters may actually be 
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a fairly gradual process. This is in contradiction to other 
authors' conclusions, which are derived from datasets with 



relati vely h i gh numbers of H^-strg ng objects (e.g. Barger 
et al. 1996; Poggianti ct al. 1999); whether this discrep 



ancy is due to selection effects or to real differences in the 
clusters studied has not yet been fully elucidated. 

A natural timescale for a gentler reduction in star forma- 
tion rate may be gleaned from the long-known observation 
that, given their present disk gas content, normal field spi- 
rals currently form stars at rates that cannot be sustained 
over a Hubble time. At face value, it appears that most 
spiral s would use up their d isk gas supply in a few Gyrs 



(e.g., Gallagher et al. 1989), though this timescale may 



be considerably extended for some choices of the IMF and 
if the effects of non-instantaneous m ass recycling is taken 
into account (Kennicutt et al. 1994). Star formation life- 
times can be significantly longer if galaxies continuously 
accre te fresh star form ation fuel from their surroundings 
(e.g., p31itz et al. 1999 ). If, as proposed by Larson, Tinsley 
& Caldwell ( |1980| ), this extended reservoir is stripped off 
a galaxy when it first enters the cluster, its star formation 
rate may decay significantly within a few Gyr, leading to 
large differences in the cluster and field populations. In 
cluste rs that are built hierarchically (as in the "bottom- 



up' scenario favored by cold dark matter cosmogonies) 
this mechanism would also establish a radial gradient in 
the star formation properties of cluster galaxies, refiecting 
the relation between the clustercentric radius of a galaxy 
and the time of its accretion into the cluster. This is an im- 
portant ingredient in the success of semianalytic models, 
which have been shown to match global cluster proper- 
ties such as the morphological composition and the blue 



1996; Kauffmann and Char- 



galax y fr action (Baugh et al. 
lot 1^ 

We explore here a simple model based on this interpreta- 
tion where the mass accretion history of a cluster obtained 
from numerical simulations of a universe dominated by 
cold dark matter is coupled with a simple star formation 
prescription for galaxies following accretion. Once cali- 
brated to reproduce the properties of local field galaxies, 
the model has no free parameters and its results can be 
compared directly with observations. We shall focus our 
analysis on a quantitative discussion of the clustercentric 
radial gradients in star formation properties and galaxy 
colors expected in this model. 

In §11 we describe the observational dataset we use, se- 
lected from the CNOCl survey. The numerical simulations 
used to derive mass accretion histories of different clusters 
are described in §|^. Our model prescriptions are presented 
in ^ and our results in We discuss the implications of 
our results in SO, and list our conclusions in SH. 



2. OBSERVATIONAL DATASET 

We use in this paper th e CNOCl cluster redshift survey 
dataset ( Yee et al. 1996| ), which consists of CFHT spec- 
tra for ^ 2000 galaxies in 15 X-ray luminous clusters at 
0.19 < z < 0.55. The observational selection effects of 
this survey are well understood a nd ar e discussed in de- 
tail in Yee, EUingson & Carlberg (1996) and Balogh et al. 
(199£). For the present analysis, we weigh the raw data 



tions which depend on the galaxy color and position on the 
CCD. One of the main advantages of this survey, especially 
for the purposes of the present study, is that the dataset 
includes spectra of foreground and background field galax- 
ies projected onto each cluster. Since both field and cluster 
galaxies are selected using the same criteria, relative dif- 
ferences between the two samples are rather insensitive to 
uncertainties in the procedure used to correct for selection 
effects. 

Since our analysis concentrates on radial gradients 
within clusters, we only consider those clusters with well 
defined centers in position and velocity space and thus ex- 
clude from the sa mple the bimodal clu sters MS0906-I-11 
and MS1358+62 ( [Carlberg et al. 1996| ). We also restrict 
the redshift range of the galaxy sample to 0.19 < z < 0.45, 
in order to facilitate comparisons with simulations ana- 
lyzed at a single epoch and to minimize effects due to 
global changes in the galaxy population as a function of 
redshift. This effectively removes two more clusters from 
the CNOCl sample, MS0016-f 16 and MS0451-03. The fi- 
nal sample has twelve clusters in total, which are scaled 
and co-added together to construct a "fiducial cluster" 
sample where effects due to substructur e and aspheric- 
ity o f individual clusters are minimized ( [Carlberg et al, 
1997[ ). The full p rocedure is described in detail in Balogh 



et al ( [1997[ ; [1999[ ). 

In brief, we u se the mass models of Carlberg, Yee & 
EUingson ( |1997 ) to divide the sample into a cluster and 
field sample; galaxies are deemed to be cluster members if 
they are within 3a of the (radially dependent) cluster ve- 
locity dispersion, and field members if they are beyond 6a. 
Cluster galaxy positions are all measured relative to the 
brightest cluster galaxy (BCG), and normalized to i?200, 
the radius at which the mean inner density is 200 times the 
critical density. The BCGs themselves are omitted from 
the final sample, as they are likely to have a unique for- 
mation history which may differ from the general cluster 
population; we briefiy compare their properties with those 
of the full sample in |5.l[ . Finally, we impose an absolute 
magnitude limit on the sample, considering only galaxies 
brighter than Mr — —18.5 -I- 51og/i at z = (Gunn-r, 
go = 0.1); when appropriately weighted, this sample is 
statistically complete. Because each individual cluster is 
at a different redshift a small evolutionary correction is 
applied to this cutoff assuming that l uminosity increase s 
in direct proportion to (1 -|- z) (e.g., Lin et al. 1999). 
At z 0.3, the luminosity cutoff we adopt is therefore 
Mr = —18.8 -I- 51og/i. This correction has little effect on 
our results because of the narrow redshift range under con- 
sideration. For the sample considered here, the luminosity 
function of the cluster population is similar to that of the 
field galaxy population. 

Individual star formation rates (SFRs) for galaxies in 
the sample are computed from the rest frame equivalent 
width of the [OIIJA3727 emission line and the rest frame 
B-band luminosity relative to solar (Lb/Lb.q) using the 
relation. 



= 3.4 X 10- 



-12/ Lb \ ^^(oiI)E(Ha)M0 yr-^ 



(1) 



by three factors: one to account for the primary selection 
effect due to source magnitude, and two secondary correc- 



Here E(Ha) is the extinction at Ha which, following Ken- 
nicutt ( |l99^ ), we take to be one magnitude. The coeffi- 
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cient in Equation ^ has been chosen so that, on average, 
the SFRs of field galaxies are consistent with their colors 
and luminosities, based on models of thei r star formation 
history that are discussed in detail in i j5.2| . This coefficient 
is consistent with t hat empirically determined by Barbaro 
& Poggianti ( 1997 ), and a bout 30% larger than that mea- 
sured by Kennicutt ( f 992 ) — well within the uncertainty 
of its determination (see Kennicutt f 99g ) . Because we are 
concerned with relative differences, neither this normaliza- 
tion nor the extinction correction has a significant effect on 
our conclusions, unless these quantities vary dramatically 
from the cluster to the field. 

Uncertainties in the equivalent widths have been as- 
sessed assuming Poisson statistics, and have been inter- 
nally calibrated to account for addit ional systematic er- 
rors, as described in Balogh et al. ( 1999| ). We exclude 
from the sample all galaxies with very large VFo(Oll) un- 
certainties, A Wo{Oll)> ISA 6% of the sample), and 
all galaxies for which [OIIJA3727 lies outside the observed 
spectral range. The final sample with measured SFRs con- 
sists of 556 cluster galaxies and 339 field galaxies. 

The W^o(Oll) measurements from which we derive SFRs 
are computed by adding up the observed flux (accounting 
for partial pixels) above the continuum level in the wave- 
length range 3713 < A/A< 3741. The continuum level is 
estimated by fitting a straight line to the flux in the ranges 
3653 < A/A< 3713 and 3741 < A/A< 3801. For weak or 
absent [Oil] features, the Wo (Oil) index will be sensitive 
to features in the continuum in these two regions. A crude 
estimate of uncertainties introduced in the Wo (Oil) mea- 
surements by these features may be obtained by using eq. 1 
to compute SFRs of galaxies expected to have little or no 
ongoing star formation; these are red cluster galaxies with 
large 4000A breaks {{g - r)o > 0.35, D4000 > 1-8) found 
within < R/R200 < 0.3. The 3cr-clipped mean SFR of 
this population is -0.057 /i~^A/0yr~^, and the standard 
deviation is 0.156 Mq yr~^. This small systematic off- 
set and uncertainty are taken into account in our modeling, 
as described in S4.2. 



3. NUMERICAL SIMULATIONS 

Cluster mass accretion rates are computed directly from 
N-body simulations of the formation of six massive clus- 
ters (0.7 < M/{1Q^^Mq) < 2.3) in a COBE-normalized, 
A = 0.7, Hq — 0.3 cosmology. The simulations are simi- 
l ar to those described in detail by Eke, Navarro, & Frenk 
We assume that "light traces mass" in the clus- 



(|T99|). 



ters and identify statistically each dark matter particle 
with a "galaxy". The virial radius, i?vir, of a cluster 
is computed using the overdensity prescription described 
in Eke, Cole & Frenk (1996), which differs slightly from 
At z = 0.3, the mean redshift of our cluster sample, 



200- 



1.4-2.4 Mpc. Each cluster has about 



R: 

Rvir ~ 1.2i?2 

9,000 particles within 2i?vir sA z = 0.3. The observations 
and model parameters are all presented in terms of the 
simulation cosmology; we use Hq — 70 km s^^ Mpc^^, 
which gives a present age of the universe of 13.5 Gyr. At 
z = 0.3 the universe is ~ 11 Gyr old. 



4. THE MODEL 



As discussed in §1, our modeling assumes that cluster 
galaxies differ from their field counterparts in their star 
forming properties because they are stripped of their sur- 
rounding gas reservoirs as they are accreted into the clus- 
ter. Star formation in cluster galaxies thus declines after 
accretion as more and more of the galaxy's remaining gas 
content gets turned into stars. 

The model involves the following steps: 

(i) All particles within 2 i?2oo from the center of each sim- 
ulated cluster at z '--^ 0.3 are selected. 

(ii) Each particle is traced back in time to find out when 
it was first accreted into the cluster. Two defini- 



tions of "accretion time" , ia 



are used. The first 



(^acc = ^cluster) is defined to be the time when the par- 
ticle ("galaxy") first finds itself within i?vir from the 
center of the current most massive progenitor of the 
cluster, and the second (tacc = igroup) is the time when 
a galaxy is first accreted into any large clump, not nec- 
essarily the most massive one. Details of this procedure 



are given in §4.1 



(iii) A SFR is chosen for each galaxy at t 



For 



simplicity, these are taken at random (using appropri- 
ate weights as discussed in §|^) from the distribution of 
SFRs computed for z 0.3 field galaxies in our sample. 
This assumes implicitly that field SFRs do not evolve 
with time, argu ably the simplest possible model. We 
explore in §5.2 the consequences of relaxing this as- 
sumption. 

(iv) We use a simple gas consumption model to estimate 
SFRs following accretion and to compute the SFR of 
each cluster galaxy at z 0.3. Details are discussed in 



4.2 



(v) Each simulated cluster is projected onto three orthog- 
onal planes. Field contamination in the ob servational 
sample is accounted for as discussed in § [l.3| . 
This procedure uniquely defines a final (z = 0.3) SFR 
for each galaxy in the simulated clusters; the mean SFR 
can then be computed as a function of clustercentric dis- 
tance, and compared with the field galaxy observations. 
We discuss now the various steps of the model in some 
detail. 

4.1. Accretion Times 

We have explored two ways of assigning "accretion 
times", iaccj to galaxies in the simulated clusters. One 
choice defines iacc = ^cluster as the time when a par- 
ticle is first found within the virial radius of the most 
massive proge nitor p resent at that time. Following Eke, 
Cole & Frenk (199q), the virial radius, i?vin is defined as 
the radius where the mean inner density of the cluster is 
p(-Rvir) = ^c{z)pc{z), where pc{z) is the critical density at 
the redshift z, and Ac(0(z)) is the "critical" overdensity 
for spherical collapse (which takes the familiar value of 178 
for = 1 models). 

The above definition of tciustor may not be completely 
appropriate, since in a hierarchically clustering universe 
particles may be first accreted into another "protocluster" 
before being accreted into the main progenitor of the final 
cluster. Therefore, we consider as an alternative definition 
of tacc the time when a particle first finds itself associated 
with a clump with circular velocity exceeding Vc = 500 
km s~^. (This circular velocity corresponds to a virial 



4 



temperature of ^ 0.8 keV). Particles are associated with 
clumps via a friends-of-friends algorithm, with an evolv- 
ing linking length parameter taken to be 10% of the mean 
inter-particle separation at each time. Accretion times de- 
fined this way are labeled tgroup- 



5.1 



that 



We shall see in 
our results are rather insensitive to the particular choice 
of accretion time definition. 

The "resolution" of iacc is limited by the number of times 
particle positions are output by the simulations. In this 
case, we have outputs every 1.34 Gyr; we therefore add a 
uniform random number between and 1.34 Gyr to each 
particle's tacc to smooth out this resolution. We do not ex- 
pect our results to be sensitive to this resolution, since we 
are combining the results of six clusters, which will smooth 
out the effect of discrete merger events. 

4.2. Star Formation Prescription 

For normal, field spiral galaxies, the SFR per unit area 
averaged over the disk, Esfr, depends on the disk gas 
content in a m anner well approxim ated b y a Schmidt law 
(Schmidt 195£). From Kennicutt ( 199§| ), this relation is 
given by 



SsFR= (2.5±0.7) X 10- 



Mopc- 



N 



Moyr ^kpc" 



(2) 



where Sgas is the average gas surface density, and the ex- 
ponent N = 1.4 ± 0.15. Converting surface densities to 
integrated values using the source diameters tabulated by 
Kennicutt, galaxies that evolve according to this relation 
and accrete no extra gas would form stars at a steadily 
decreasing rate given by, 



model SFR. These parameters were chosen from the SFR 
distribution of red, central cluster galaxies (with h — 0.7), 
as described in §^ 

4.3. Field Contamination 

Observational cluster datasets are contaminated by field 
galaxies, projected onto the cluster, that have Hubble-flow 
redshifts similar to the peculiar velocities induced by the 
cluster potential. Our model takes this into a ccount follow- 
ing the procedure of Carlberg et al. (1997). In brief, the 
velocity differences between cluster and individual galax- 
ies in each dataset are scaled to the velocity dispersion 
of each cluster and their clustercentric projected radii to 
i?200- These normalized velocities, VJiorm — Ati/cr, and 
radii, Rpmj/ Rioo^ are then combined into a single dataset. 
The density of galaxies per unit Korm at 5 < Korm < 25 is 
computed in radial bins to account for the radially vary- 
ing (7. Under the assumption that this density remains 
constant within each bin, it is used to calculate the ex- 
pected number of interloper galaxies, i.e., those within 3<t 
in each radial bin. This provides a direct estimate of the 
fraction of galaxies deemed cluster members that are ac- 
tually field galaxies projected onto the cluster in "redshift 
space". This fraction is typically ^ 1% in the central clus- 
ter regions, but may be as large as ^ 12% at i?200- Our 
modeling accounts for this effect by including an appro- 
priate number of "field" galaxies (randomly selected from 
the field sample) in all computations. 



5. RESULTS 



5.1. 



M4t') =M*(0) ( 1 + 0.33- 

tr 



-3.5 



M, 



(3) 



Non- evolving field SFR model 

Figure |l| shows the mean SFR per galaxy as a function of 
projected distance from the cluster center, normalized to 
i?,200- The solid squares correspond to the CNOCl data for 
the complete sample of galaxies brighter than Mr = —19.6 
(at z = 0.3, h = 0.7), averaged over radial bins and with 
l-CT jackknife error bars. The mean SFR per galaxy in- 
creases systematically from the center of the cluster out- 
wards, from almost zero near the center to about ~ 0.7 Mq 
yr~^ in the outskirts of the cluster. We note that the ab- 
solute SFR values quoted here are sensitive to a number 
of sample selection parameters; in particular to the lumi- 
nosity cutoff and the uncertain coefhcient in eq. ^ so care 
must be exercised when comparing these results to other 
work. On the other hand, relative differences between clus- 
ter and field populations are robust, since the two samples 
have similar luminosity functions, and are drawn from the 
same survey with identical selection criteria. 

As shown in Figure I, the average SFR per galaxy in 
the field is significantly higher than in clusters. Remark- 
ably, even at radii as far from the center of the cluster 
as 2i?2oo, cluster star formation rates remain depressed 
by almost a factor of two relative to the field.*. We note 
here that the BCGs, which are omitted from the data sam- 
ple, have unusually high SFRs; the mean SFR of the 11 
BCGs which satisfy our selection criteria^ is 14 AIq yr~^, 

^This large factor of R approximates the effects of a non-instantaneous recycling calculation, assuming a Scale IMF, and is fairly good for 
slowly evolving disks; for rapidly evolving disks, _R can be as high as 0.75. 

^Because the luminosity functions of the field and cluster population considered in this sample are quite similar, the same result is obtained 
for the SFR per unit luminosity as for the SFR per galaxy. 

5No spectrum is available for the BCG in MS 0451.5-1-0250. 



where (0) is the initial SFR and te ~ 
1.48 (Af*(O)/Af0yr"i)~"-29 Gyr is the characteristic gas 
consumption timescale. As discussed in we adopt 
M*(0) values taken at random from the measured field 
SFRs and set t' = t — tacc- With these assumptions it is 
possible to compute SFRs of cluster galaxies at z 0.3, 
the mean redshift of the observational sample. In the 
above calculation, we have neglected the effects of gas 
recycling. If we assume instantaneous recycling, where a 
fraction R of every solar mass of stars formed is returned 
as gas, the timescale tp is increased by a fact or 1/(1 — R). 
Assuming R = 0.33 ( Kennicutt et al. 1994 )*^, we obtain 
te « 2.2(M*(O)/M0yr-i)-O-29 Gyr. This increase does 
not have a significant effect on the results we present, so 
we neglect the effects of recycling throughout this work. 

Finally, we correct the values obtained from eq. ^ for 
the systematic offset and uncertainty in the SFR measure- 
ments discussed in This is done by adding a random 
number, drawn from a Gaussian distribution with a mean 
of — O.I2M0 yr~^ and a variance of 0.32M(?, yr~^, to the 
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much greater than the mean field value. Only four of these 
eleven galaxies have no significant [Oil] emission. The low 
redshift BCGs are discussed in more detail in Davidge & 
Grinder ( [1995D ; we do not consider these galaxies further 
in this study. 

Figure ^ also shows, with open symbols, the results of 
the modeling procedure outlined in f^. Open squares and 
triangles correspond to the two different accretion time 
definitions discussed in §4.1. Error bars represent the l-tr 
variance of the 18 numerical realizations (6 simulated clus- 
ters and 3 orthogonal projections per cluster). The agree- 
ment between the model and the observations is remark- 
able, especially considering that the modeling involves no 
free parameters. The model reproduces the observed SFR 
gradient and even the observed depression, relative to the 
field, of SFRs outside i?2oo • 

The latter result is somewhat surprising, since in spher- 
ical accretion models particles outside R2m are infalling 
into t he clu ster for the first time ( Bertschinger 1985; White 
et al. 1993) and therefore their SFRs have yet to feel the 
effects of the cluster environment. The main reason for our 
result is that a substantial fraction (54 ± 20% for the six 
clusters we studied) of particles between 1 and 2 i?2oo have 
actually been inside the virial radius of the main progen- 
itor at some earlier time. These are often particles that 
populated the outskirts of recently accreted clumps and 
that, although still bound to the system, have been scat- 
tered to large apocenter orbits during the merger process. 

Interestingly, assuming that the onset of the SFR de- 
cline occurs when a galaxy is accreted into any clump 
with circular velocity exceeding 500 km rather than 
the cluster's main progenitor has only a small effect on 
this result within i?200 (witness the good agreement be- 
tween open triangles and squares in Figure ||). Note that 
the mean SFR beyond i?200 is further suppressed under 
this assumption, resulting in even better agreement with 
observations. This is because some particles beyond i?200 
are found within fairly massive groups, although they may 
never have been within the virial radius of the main clus- 
ter. 

An important feature of this model is that many clus- 
ter galaxies at z = 0.3 have substantial SFRs. In par- 
ticular, near R200: ^ 20% of cluster galaxies have SFRs 
in excess of IMQyr"^, and this declines to about 10% 
at i? = 0.5i?2oo; fractions which are consistent with the 
CNOCl data. These large SFRs are not the result of 
cluster-induced starbursts, but correspond to recently ac- 
creted field galaxies in which star formation has not yet 
been co mplet ely quenched. In a recent study, Balogh & 
Morris (200C) have measured Ha equivalent widths for 
galaxies in Abell 2390, and they failed to find a substan- 
tial population with large Ha fluxes that were undetected 
in [Oil]. Thus, there does not appear to be a population 
of dust-obscured starburst galaxies in this cluster which 
were missed in the CNOCl survey. Our modeling indi- 
cates that cluster-induced starbursts are not necessary to 
generate the levels of star formation seen in these clusters. 

We conclude from this comparison that a model based 
on the simple assumption that the SFR of a galaxy de- 
creases on gas consumption timescales is able to reproduce 



the data extremely well, lending support to the underlying 
hypothesis that continuous accretion of external gas is re- 
sponsible for maintaining the SFRs of normal spirals over 
a Hubble time. 



5.2. Evolving field SFR model 

There are two weaknesses in the model explored in the 
previous sections. One is the assumption that the SFRs of 
field galaxies, which are used to assign "initial" SFRs to 
cluster galaxies at accretion time (see eq. p|) , do not evolve 
with time. Given the mounting evidence that the average 
SFR per unit volume evolves strongly w ith lookback time 
( [Lilly et al. 1995| ; [Madau et al. 1996] ) it is necessary to 
explore the consequences of relaxing the non-evolving SFR 
assumption on our results. The second is the luminosity 
evolution of cluster galaxies, which must also be modeled 
in order to account for galaxies that may fade beyond the 
observational magnitude limit when their SFR declines. 

A proper treatment of these effects, which must take into 
account the merger and star formation history of galax- 
ies, is of great interest, but well beyond the scope of this 
study. In order to at least explore the effects that this more 
complete treatment will have, we use simple t— models to 
model the SFR evolution of field galaxies in our sample, 
i.e., by assuming that, 



Ai(t) = ^i(^o)e-(*-*°)/^ 



(4) 



where M,(to) corresponds to the SFR of galaxies in the 
field at t = ip ~ 11 Gyr, the age of the universe at z ~ 0.3. 
Once T is determined for each galaxy it is possible to con- 
struct, at arbitrary times, an "evolving" field SFR distri- 
bution which matches the observations at z ~ 0.3. For 
a given IMF, and constant reddening, r is uniquely de- 
termined by the observed {g — r)o color ^ of the ga laxy, 
as shown in Figure 0. We ha ve used the PEGASE (Fioc 



and Rocca-Volmerange 1997) spectrophotometric models 



with a Salpeter IMF (with lower and upper limit given by 
0.1 < M/Mq < 120) to determine r for each galaxy in 
our field sample. The bluest galaxies {{g — r)o < 0.04) 
have T < 0, corresponding to a SFR that increases with 
time, while moderately blue galaxies {{g — r)o ~ 0) have 
T approaching (positive or negative) infinity, correspond- 
ing to a constant SFR. The reddest field galaxies, those 
with r < 300 Myr, are assumed to have formed in a single 
burst at high redshift, and are modeled as such in the clus- 
ter. The net result is a rather extreme model where the 
global star formation rate in the field population increases 
steadily out to z w 10. We explore below as well the conse- 
quences of restricting this increase at large z, and conclude 
that our results are quite insensitive to the precise nature 
of the redshift evolution in the field. 

The normalization constant of cq. ^, Af,(to), may in 
principle be determined by the observed field galaxy SFR 
at but we found that, because of the large uncertainties 
in individual SFR determinations, this procedure leads to 
large discrepancies in the total luminosity of galaxies in the 
sample once the SFRs are integrated over time. Because 
total luminosities are considerably better determined than 



®This color is available for all galaxies in the CNOCl dataset. We have made a small adjustment to the zero point of the observed [g — r)o 
colors, reducing them by 0.05 so that the reddest observed galaxies have the colors predicted by a model in which all stars formed at t = in 
an instantaneous burst. 
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Fig. 1. — The mean SFR per galaxy as a function of projected radius for galaxies in the CNOCl cluster sample 
(solid squares) compared with the model predictions under the assumption that iacc = ^cluster (open squares) and that 
iacc = igroup (opcn triangles). The horizontal dashed lines correspond to the field SFR, bracketed by its 1-cr dispersion. 
The SFR gradient in the CNOCl clusters is accurately reproduced by our simple accretion models. Error bars are all la. 
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Fig. 2. — The histogram shows the distribution of (g — r)o colors for the field galaxy sample, in arbitrary units. The solid 

rmbols correspond to model colors at z = 0.3 (an age of 11 Gyr) for galaxies with star formation history given by eqn. 
and parameter t as shown on the ordinate axis. The solid line through these points is the linear interpolation used to 
assign r to each galaxy, based on its observed {g — r)o- 
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SFRs in our dataset, and are less sensitive to burst-like 
events that eq. ^ cannot reproduce, we compute M^,{tQ) 
for each field galaxy by requiring the model to reproduce 
its observed total Gunn-r luminosity at z — 0.3. This 
choice implies that, on a galaxy by galaxy basis, the model 
z = 0.3 SFRs no longer correspond to those derived from 
the M^o(OII) measurements. The observed and model field 
SFR distributions are compared in the bottom right panel 
of Figure ||; where the coefficient in eq. 
has been chosen so that both distributions have identical 
averages. Although the difference is statistically signifi- 
cant, the two distributions are qualitatively fairly similar, 
indicating that, on average, the simple r-model provides 
a reasonable description of the z = 0.3 field. This con- 
sistency allows us to compare model and observed SFRs, 
even though the former is based on galaxy color, while the 
latter is determined from nebular emission. 

A further advantage of modeling the SFR history of the 
field population is that it allows us to model the color evo- 
lution of cluster galaxies. Color {{g — r)o) distributions of 
model galaxies (which include the addition of a Gaussian 
random number to account for the 0.02 magnitude error 
typical of the CNOCl data) are compared with the obser- 
vational data in the left panels of Figure H The model 



shown adopts ta 



^group, and excludes galaxies fainter 



than Mr = —19.6 aX z = 0.3. This model is very success- 
ful at producing a cluster core dominated by red galaxies 
from an initial field galaxy sample that has a much broader 
and bluer color distribution. In the outskirts of the clus- 
ter, both the observed and model color distributions show 
a significant increase in the population of blue galaxies. 
The qualitative agreement between model and observed 
SFR and color distributions throughout the cluster is quite 
remarkable for such a simple prescription. 

The right panels of Figure ^ compare the observed SFR 
distributions with this same model. In agreement with ob- 
servations, model distributions show a strong reduction in 
the fraction of galaxies forming stars at rates higher than 
about ~ IM0 yr~^ near the cluster center. The agreement 
between model and observations is quite good, although 
model galaxies form stars at slightly higher rates in both 
bins. 

The differences between model and observed mean clus- 
ter SFRs are clearly apparent in Figure ^, which shows the 
clustercentric gradient predicted by three variants of the 
evolving field-SFR model. The open circles correspond to 
the fiducial model as described above, and are generally 
larger than the observed cluster mean. The open triangles 
show how this result changes if we neglect galaxy fading; 
i.e., when we include all cluster galaxies in the average. 
This recovers the good match to the data, and shows that 
it is primarily the luminosity fading of galaxies that affects 
our result, and not the nature of the redshift evolution it- 
self. To demonstrate this explicitly, we present a model 
(open squares) in which the field galaxy SFRs are evolved 
as in the fiducial model but are held constant for z > 1.5. 
The results obtained with this model do not differ signif- 
icantly from those of the fiducial model, lending support 
to our conclusion that our results depend only mildly on 
the SFR redshift evolution. 

We wish to stress that the purpose of the above exer- 
cise is to assess the sensitivity of our modeling to various 



assumptions about the redshift dependence of the SFR of 
the field galaxy population rather than to build a real- 
istic model of field galaxy evolution. We conclude that 
the luminosity evolution of the field galaxies has a small, 
but non-negligible effect on our results. Although it ap- 
pears that qualitatively our conclusions are safe, it is clear 
that definitive answers will have to wait for a more realis- 
tic modeling of the field galaxy star formation evolution, 
such as the one imple mented in semiana lytic models of 
galaxy formation (e.g. Baugh et al. 1996, Diaferio et al. 
in preparation). 

6. DISCUSSION 

We present models of the clustercentric dependence of 
star formation and colors of galaxies in rich clusters, under 
the following assumptions: (i) the cluster galaxy popula- 
tion is built by the ongoing accretion of field galaxies, (ii) 
SFRs in field galaxies are sustained by regular accretion 
of gas from their surroundings, and (iii) reservoirs of fresh 
star formation fuel are lost as galaxies plunge into the clus- 
ter potential. 

Within this context, our model provides support for a 
gradual decline (over a timescale of a few Gyrs) in the star 
formation rates of cluster galaxies after accretion. Actu- 
ally, results similar to those presented in the previous sec- 
tion may be obtained if the SFR in all cluster galaxies is 
assumed to decay exponentially after accretion with fixed 
timescales 1 <J 3 Gyrs. Decline timescales longer than 
^ 3 Gyr lead to unacceptably large star formation rates 
in model clusters at z ^ 0.3. On the other hand, sharp 
truncation of star formation (tg <^ 1 Gyr) would result in 
too little star formation within clusters and, furthermore, 
would lead to an abundant population of galaxies with 
strong Balmer lines but no nebular emission lines (K-l-A 
galaxies), and these appear to be r are in the very lum i- 
nous X-ray clusters we study here (Balogh et al. 1999). 
Larger fractions of K-fA galaxies have been reported in 
other cluste r datasets, in particu l ar by the MORPHS co l- 
laboration (Dressier et al. 199£; Poggianti et al. 1999), 



but it is still unclear whether this apparent disagreement 
is a result of the procedure used to select the spectro- 
scopic sample, the effects of dust obscuration, or perhaps 
a genuine effect of the dependence of SFRs on other clus- 
ter properties such as X-ray luminosity or temperature 
( [Balogh et al. 1999 ). Further analysis is required in or- 
der to assess whether the simple model we propose here 
is applicable to clusters other than the relatively massive. 
X-ray luminous systems targeted by the CNOCl survey. 

With this caveat, the success of our model strongly sug- 
gests that (i) gradients in galaxy properties arise from 
gradients in accretion times and (ii) that the cessation of 
star formation need not take place abruptly to explain the 
observed SFRs of cluster galaxies. On the other hand, 
our understanding of the star formation history of cluster 
galaxies is bound to remain incomplete until the physical 
mechanism responsible for the decline of star formation in 
cluster galaxies is fully elucidated. The loss of external gas 
reservoirs advocated here is attractive, but conclusive ob- 
servational evidence that such reservoirs exist in isolated 



field gal axies has been slow to emerge (Benson ct al. 1999 
but see Blitz et al. 1999). Other processes that may re- 



duce star formation rates, such as tides, harassment, and 
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Fig. 3. — Left panels: {g — r)o distributions in the observed CNOCl sample {heavy solid lines) and our evolving field-SFR 
model (thin solid lines). The samples arc divided into a field sample and two cluster samples (inner and outer regions), as 
labeled in each panel. By construction, the model color distribution matches the field observations exactly. Right panels: 
SFR distributions for the same samples shown in the left panels. The model provides a good match to the color and SFR 
distributions of the observed cluster population, although it slightly overestimates the mean cluster SFR. 
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Fig. 4. — Same as Figure |i| but for three evolving field SFR models. In the fiducial model (open circles), the SFRs in field 
galaxies evolve as dictated by simple t models normalize d to match the colors of field galaxies at z ~ 0.3, and galaxies 



that fade below the luminosity limit are excluded (see §5^). Accretion times are defined as iacc = tgroup. The open 
squares show the result of modifying the evolution such that each galaxy's SFR remains constant at z > 1.5. Neglecting 
the effects of galaxy fading in the fiducial model results in the open triangles. The data are shown by the filled squares, 
connected by the solid line. 



ram pressure tripping, operate on similar timescales once 
a galaxy has been accreted into the cluster and it is there- 
fore unlikely that analysis of the kind we present here will 
be able to distinguish clearly amongst them. 

Another major question that our models do not address 
is the origin of gradients in galaxy morphology that par- 
allel the color and s tar formation gradients we focus on 
(e.g., Dressier 1980 ; Dressier et al. 1997 ; Koopmann ancj 
Kenney 199S). The over-representation of ellipticals near 
the center of X-ray luminous clusters may reflect higher 
merger rates between nearly equal mass systems in sys- 
tems that collapse early to form cluster cores, but this is 
an issue that remains unexplored in our model. The con- 
struction of large spectroscopical datasets that probe the 
dependence of galaxy SFRs, colors and morphologies on 
cluster properties such as concentration, richness, velocity 
dispersion and X-ray properties, coupled with numerical 
and/or semianalytical models that treat self-consistently 
the accretion history and dynamics (and hydrodynamics) 
of galaxy and cluster formation, are probably the most 
promising way to make substantial progress in the sub- 
ject. 

7. CONCLUSIONS 

We present a simple model to account for the systematic 
differences in the star formation properties of galaxies in 
CNOCl clusters and the field. The model assumes that 
the cluster is built through the accretion of field galaxies 
whose star formation rates gradually decline after entering 
the cluster as a result of the removal, through tides or ram 
pressure, of the gaseous envelopes needed to supply nor- 



mal spirals with the fuel needed to sustain their present 
star formation rates over a Hubble time. Once calibrated 
to reproduce observations of nearby spirals the model has 
no free parameters. 

Using cluster mass accretion rates determined from N- 
body simulations of cluster formation in a ACDM universe, 
our model is able to reproduce qualitative and quantita- 
tive differences in the mean star formation rates and colors 
between clusters and the field. The model demonstrates 
that the origin of radial gradients in these properties is 
the natural consequence of the strong correlation between 
radius and accretion times which results from the hierar- 
chical assembly of the cluster. Interestingly, the model also 
explains why star formation in the outskirts of clusters is 
found to be almost a factor of two below the field average 
as far out as twice the virial radius of the cluster. This is 
a result of cluster members being pushed onto highly ec- 
centric, loosely bound orbits during major merger events. 
Our results are robust to evolution in the star formation 
properties of field galaxies, but somewhat more sensitive to 
the manner in which fading of cluster galaxies is modeled. 
Realistic modeling of the field population star formation 
histories is required to fully understand the consequences 
of this fading. 

We conclude that the stripping of extended gaseous 
reservoirs by the cluster environment and the gradual de- 
cline that follows gas consumption is likely to be the main 
mechanism that differentiates the star formation proper- 
ties of cluster and field galaxies. 

We have made extensive use of the CNOCl dataset 
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